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BEYOND MEAN FIELD: ON THE ROLE OF PAIR 
EXCITATIONS IN THE EVOLUTION OF 
CONDENSATES. 

M. GRILLAKIS AND M. MACHEDON 

This paper is dedicated to Professor Choquet-Bruhat, whose work on nonlinear 
wave equations was an inspiration to us earlier in our careers 


1. Introduction 

Abstract. This paper is in part a summary of our earlier work 
mmm, and in part an announcement introducing a refinement 
of the equations for the pair excitation function used in our previ¬ 
ous work with D. Margetis. The new equations are Euler-Lagrange 
equations, and the solutions conserve energy and the number of 
particles. 


2. Introduction 


The problem, which has received a lot of attention in recent years, is 
concerned with the evolution of the iV-body linear Schrodinger equation 


ld_ 
i dt 


ip N (t , •) = H N 'ip N (t , •) with 


f/hv(0, xi, • • ■ , x N ) = (p 0 (x 1 )(p 0 (x 2 ) ■ ■ ■ 4>o(x n ) 


||') ||l 2 (R 3JV ) = 1 


The Hamiltonian is an operator of the form 


Hn = Y1 Ax > ~ VN ( Xi ~ 

3 =1 i<j 

where v^{x) := N 3/3 v(N l3 x) with 0 < ff < 1 models the strength of two 
body interactions. Notice that if f3 > 0 then vn(x) —> 8(x) as N — > oo. 


The authors thank Sergiu Klainerman for the interest shown for this work, and 
John Millson for many discussions related to the symplectic group and its represen¬ 
tations. The authors would like to thank the Department of Applied Mathematics 
at the University of Crete and ACMAC for their hospitality during the preparation 
of the present work. 
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For simplicity we assume that v E Co and v > 0. The goal is to show, 
in a sense to be made precise, 

xi, ■ ■ • , x N ) ~ e lNx{t) (f)(t, xi ) 0 (t, x 2 ) ■ ■ ■ 0 (t, x N ) ( 1 ) 

where 0 satishes a suitable non-linear Schrodinger equation. In partic¬ 
ular, this approximation is not true in L 2 (R 3Ar ). 

The motivation for this problem is that in the presence of a trap the 
ground state of H N looks like 

*iv(Zl, X 2 , ■ ■ ■ , X N ) ~ 0 O (Xi)0o(x 2 ) • • • 00 (x N ) 

This is suggested by the result of Lieb and Seiringer who showed in 
[ 2 B] that 

li{ x i x> ) -t M X )M X ') 

where 

qf ( x , x ) = I ^ n( x , x 2j , X N )ty N (x', x 2i ' ' ' j x N )dx 2 ■ ■ ■ dx N 

Here ||0o||l 2 = 1 and 0o minimizes the Gross-Pitaevskii functional. See 
[25] for extensive background. 

The reason for the recent attention to this problem is two-fold. On 
the one hand experimental advances during the last twenty years made 
the creation and manipulation of condensates in the laboratory possi¬ 
ble, on the other hand recent mathematical developments made possi¬ 
ble the rigorous treatment of the equations when the number of parti¬ 
cles, namely N, is large. 

While this is a ’’classical PDE problem” (as opposed to a Fock space 
problem), the PDE approach to this problem has only been studied 
systematically during the last 10-15 years, in the series of papers of 
Erdos and Yau [8j, and Erdos, Schlein and Yau [9] to [Uj. See also [7j. 
These papers prove 

7 i(t, cc, x> ) ->• 0 (t, x) 0 (t, x') ( 2 ) 

in trace norm as N —» oo, and similarly for the higher order marginal 
density matrices 7 ^, where k is fixed. The problem becomes more 
difficult and interesting as the parameter /3 in the definition of Vn 
approaches 1. The strategy of these papers is based on the older work 
of Spohn [30]. Recent simplifications and generalizations, based on 
harmonic analysis techniques and a ’’boardgame argument” inspired 
by the Feynman diagram approach of Erdos, Schlein and Yau, were 
given in [ 21 ], [ 22 ], [ 6 ], [3], [I], [5]. See also [H], [27] for a different 
approach. 
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The symmetric Fock space approach to the problem is much older. 

It originated in physics, with the papers by Lee, Huang and Yang [23] 
in the static case, and Wu [31] in the time-dependent case. See also 
[2]. It continued with the mathematically rigorous work of Hepp [20] . 
and Ginibre and Velo [15]. 

Motivated by the goal of obtaining a convergence rate to solutions of 
NLS in ([2]), Rodnianski and Schlein resumed the rigorous Fock space 
approach in [28] • This paper, as well as the older work of Wu, served 
as an inspiration for our work. Our goal is to obtain a refinement to 
Q which provides an L 2 (M 3JV ) and Fock space estimate. This leads to 
the introduction of the pair excitation function k. 

We also mention the recent preprint [1] where a similar approach 
(but with an explicit choice of pair excitation function k ) is used to 
prove convergence of the density matrices in the critical case (3 = 1. 

3. Fock space 

In this section we briefly review symmetric Fock space, following the 
notation of [!§]• See [28], for more details. The elements of F are 
vectors of the form 

\^) = {4’0 , -01 Og) , i>2{X\,Xi) , ... ) 

where -0o £ C and 'i/’fc are symmetric L 2 functions. The norm of such a 
vector is, 

OO 

II WWt = W 7 ^) = l^°l 2 + IWII2 • 

n =1 

The creation and anihilation distribution valued operators denoted by 
a* and a x respectively which act on vectors of the form (0, • • • , 0 n _ 1, 0, • • •) 
and (0, • • • ,0 n+ i,O, • • •) by 

1 n 

a*(-0 n _i) := — 7 = ^ 5(x - Xj)-0 n _ i(xi,..., xj- 1, x j+1 , ...,x n ) 

V ^ -1 

v J= i 

«x(0n+1) := VnTli> n+ i{[x\,x 1, ...,x n ) 

with [x] indicating that the variable x is frozen. The vacuum state is 
defined as follows: 

| 0 ) := ( 1 , 0 , 0 ...) 

and aa,|0) = 0. One can easily check that [a x , a *] = S(x — y) and since 
the creation and anihilation operators are distribution valued we can 
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form operators that act on F by introducing a field, say 0(x), and form 


a(0) := 


dx {cj)(x)a x } 


and a* (0) 


J dx {4>{x)a* x } 


where by convention we associate a with 0 and a* with 0. These 
operators are well defined, unbounded, on F provided that 0 is square 
integrable. The creation and anihilation operators provide a way to 
introduce coherent states in F in the following manner, first define the 
skew-Hermitian operator 

*4.(0) := J dx {(j)(x)a x — 0(x)a* j (3) 

and then introduce iV-particle coherent states as 

|0(0)) : =e -VAU(*)| 0 ) . ( 4 ) 

This is the Weyl operator used by Rodnianski and Schlein in [2Ej. It 
is easy to check that 

e -VNA(4>) | 0 ) = ^ Cnfl^Xj) ... j with c n = (e~ N N n /7i\) 1/2 . 

In particular, by Stirling’s formula, the main term that we are inter¬ 
ested in has the coefficient 


C/v ~ (27T N) X / 4 


(5) 


Thus a coherent state introduces a tensor product in each sector of F. 

For the construction analogous to (J3]) involving quadratics, start with 
the Lie algebra of real or complex symplectic ’’matrices” of the form 

T _ (d(x, y) l(x,y) \ 

L -{k(x,y) —d T (x, y) J 


where d, k and l are kernels in L 2 , and k and l are symmetric in (x, y ). 
We denote this Lie algebra sp( C) or sp(M) depending on whether the 
kernels d, k and l are complex or real. The natural setting for us (which 
will insure that the Fock space operator e 1 ^ defined below, is unitary, 
see also the appendix of [18]) is the subalgebra sp c (M) = Wsp(M)W -1 
where 


W = 


1 

7 ! 


The elements of sp c (R) look like 




fid(x,y) k(x, y) \ 
\k(x,y) —id T (x, y)J 


(6) 


5 


with L 2 kernels d complex and self-adjoint, k complex and symmetric. 


Remark 3.1. The corresponding group elements E G Sp c (R) (in par¬ 
ticular E = e L , L G sp c (M)) satisfy the following three properties: 

• E commutes with the real structure a defined by cr(0, 0) = 
(0, 0), in other words E is of the form 

F= (P( x ,y) Q(x,y)\ 

\Q(x,y) P{x,y)J 


• E belongs to the infinite dimensional analogue of U(n, n) , in 
other words 


E* 






• E is in the symplectic group, meaning 


E t 






In fact, any two of the above imply the third. The conceptual rea¬ 
son for this is that the symplectic inner product ^(0i,0i), (02,02= 

/ 0102 — / 0102 and the ”C/(n, n)” inner product ( ( 0 i, 0 i), ( 0 2 ,02) ) = 


/ 0102-/0i02 are related by ^( 0 i, 0 i), (02, 0 2 )y = j^( 0 i, 0 i), a( 0 2 , 0 2 ) 

See Folland’s book ra for more along these lines in the finite dimen¬ 
sional case.These matrices are called Bogoliubov rotations in [1], 


Our approach is based on the map from L G sp( C) to quadratic 
polynomials in (a, a*) in the following manner, 


m = i 




(7) 


- J dxdy {d(x,y)a x a* + d(y,x)a*a y + k(x,y)a*a* — l(x,y)a x a y } 


This is the infinite dimensional Segal-Shale-Wcil infinitesimal repre¬ 
sentation. The group representation was studied in [2pj. The crucial 
property of this map is the Lie algebra isomorphism 

[P(Li),I(L 2 )\ = X{[L 1 ,L 2 ]) (8) 

Notice that if L G sp c (M), then L has the form (0 and Z(L) is skew- 
Hermitian, thus eXO is a unitary operator on Fock space. For the 
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applications that follow we will only use the self-adjoint elements of 

spJM) 


K = 


0 k(t, x, y) 
k(t,x,y) 0 


(9) 


and the corresponding 


m :=X(K) = - 


j dxdy { k(t, x, y)a x a y — k(t, x, y)a*a*} . (10) 


( 0 k(t,x,y)\ 

\k{t,x,y) 0 ) 


We easily compute 


k = ( ch(fc) sh (k) \ 

\sh(/c) ch(/c) / 

where 

sh (k) := k + —k o k o k + ... , 
ch(k) := S(x - y) + ^ k o k + ... , 


( 11 ) 


(12a) 

(12b) 


This particular construction and the corresponding unitary operator 
e B were introduced in [17]. 

The Fock Hamiltonian is 


'H\='H 1 - N~ l V where, (13a) 

Hi := j dxdy {A x S(x — y)a*a y } and (13b) 

V := ^ J dxdy {uat(x - y)a* x a* y a x a x } . (13c) 

It is a diagonal operator on Fock space, and it acts as a regular PDE 
Hamiltonian in n variable 

n 

H„, PDE = E N, - ^ E - **)) 

3 —1 


on the nth component of F. 
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4. Outline of older results 


Our goal is to study the evolution of coherent initial conditions of 
the form 

|t/w) = e ^e-^^ o) |0> (14) 

The papers [17, 181 HHj propose an approximation of the form 

(15) 

and derive Schrodinger type equations for (p{t,x), k(t,x,y ) so that 
exact(t)) ~ e lNx ^\'ipapprit)), with y(f) a real phase factor, and find 
precise estimates in Fock space, see Theorem (14. ip below. Our strategy 
is to consider 

\Aed) = e B We^)e^e-^(o)| 0 > 

and then find a ’’reduced Hamiltonian” H re d so that 

~dt\^red) T'i'red | ^red) • ( 16 ) 

7/ 

The reduced Hamiltonian is 

Hred ■= -(< d t e B )e~ B 

i v 7 

+ e B (^{d t e^ A )e~^ A + e^ A Ue~^ e" 6 . 

It can be written abstractly as a composition (in space only) of opera¬ 
tors 

nj - 1 ^ I ( 1 ^ , nj\ r, — y/NA^—B 

Hr “‘-lm +ee {~lm +n ) oe e 

Explicitly it is 

Hed = NV 0 + N l ' 2 e B V ie ~ B 

+ He + T{R) - N~ 1 / 2 e B V 3 e~ B - fWTV 5 ( 17 ) 

where the various terms are defined below. V n indicate polynomials of 
degree n in a, a* to be given explicitly: 



~7;J dxdy {v N (x - y)\(j)(x)\ 2 \(j)(y)\ 2 } . 

+ K} 

=a*(h(t , ■) + a(/i(l, ■)) 


(18) 

( 19 ) 
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where h := —(1 /i)d t (j) + A0 — (v^ * |0| 2 )0. 


7-Lg : = - J dxdy {-g N (t,x,y)a* y a x - g N (t,y,x)a* x a y ] (20a) 

where (20b) 

g N (t, x , y) := -A J(x - y) + (v N * |0| 2 )(t, x)5(x - y ) 

+ ujv(ic-j/)0(t,a;)0(i,2/) (20c) 


and 


i? = - (^-e K \ e~ K + [G, e K ]e~ K + e A 'Me.- A ' = 
i \ot J 

_ f-W(ch(k)) -S(sh (k))\ o / ch(fc) f21 , 

\ S(sh(fc)) W(ch(fc)) J \-sh(k ) ch(fc) y 

+ e K Me~ K 


where S describes a Schrodinger type evolution, while W is a Wigner 
type operator by 


S (s) := 
while 
m(x, y ) 


-s t + g T o s + s o g and W(p) := -p t + [g T ,p\ 
X X 

M( ° where 

V — m 0 I 

:= —v N (x — y)(ft(x)(J)(y), v N (x) = N 3/3 v(N l3 x) 


and G:=(/ _° T ) . 


Finally, 


V, := 


[A, V] = I dxdy {v N (x - y) ((p(y)a* x a* y a x + cj>(y)a* x a x a y ) } 


V A := V = (1/2) / dxdy {v N (x - y)a* x a* y a x a y } . 


(22a) 

(22b) 


The main result of (19] , building on the previous papers of the authors 
and D. Mar get is HUES], can be summarized as follows. 
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Theorem 4.1. Let 0 and k satisfy 

-d t 4> — A0 + ( 'v N * |0| 2 )0 = 0 (23a) 

and either one of the following equivalent equations: (23b) 

1) (sh (k)) — ch(/c) o o ch(fc) = ^ch(/c)j + sh(fc) orn^j o sh (k) 


(23c) 

or else the equivalent non-liner equation (23d) 

2) S(th(fc)) = m + th(A;) onto th(fc) (23e) 

where th(fc) := ch (k) o sh (k) (23f) 

or else the equivalent system of liner equations (23g) 

за) S (sh(2/c)) = m,N ° ch(2fc) + ch(2fc) o (23h) 

зб) W ^ch(2/c) j = rn jv o sh(2 k) — sh(2/c) o Win • (23i) 

with prescribed initial conditions 0(0,-) = 0o, /c(0, •,•) = 0. If 0, k 


satisfy the above equations, then there exists a real phase function x 
such that 

III- e w *<'>|^(()>|U < C(1+ A ,?,!.w ( 2 1 + t) ■ (24) 

provided 0 < 0 < |. 

The purpose of the present paper is to introduce and study a coupled 
refinement of the system fj23aH . (I23hj) . (I23il) which, we believe, is the 
correct system describing the case 0 = 1. These equations occur as 
Euler-Lagrange equations, and are written down explicitly in Theorem 

dSU). 


5. Main new results 

Since Hred is a fourth order polynomial in a and a*, 

Hred\0) = (X o ,X 1 ,X 2 ,X 3 ,X 4 ,0,---). (25) 

Definition 5.1. Define the Lagrangian 

£ = — J X 0 (t)dt (26) 

The new, coupled equations for 0 and k that we introduce in this 
paper are A0 = 0 and X 2 = 0. 
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We first prove that C is indeed the Lagrangian for these equations. 
We start by showing ’’abstractly” that 


= v'N J ch(k)(t,x,-) 

^ = -Ldlp) ° W o ch(fc) 


Xi(t, x)sh(k)(t, x, y)) dx (27) 

(28) 


where ( = th(fc) = ch(fc) o sh (k). 
zeroth order term X 0 (t) in l-i re d |0) 
density for our coupled equations): 


We then compute explicitly the 
(which provides the Lagrangian 


~X 0 (t) 




dx 1 dx 2 Vi_ 2 \(j) 1 (j )2 + — (sh o ch)i )2 
dxidx 2 dx 3 v^_ 2 \(j) 1 sh2,3 + </> 2 sh 1)3 | 2 


2 


+ 


1 

2 





(shi^^shi^) + (Vi^shi^ 


1 

2N 


dtdxidx 2 Vi_ 2 l | (sh o sh)i j2 | 2 + (sh o shji^sh o sh) 2]2 


where shi, 2 is an abbreviation for sh(k)(t, x\, x 2 ), v^_ 2 = Vn(x\ — x 2 ), 
etc, and the products are pointwise products, while compositions are 
denoted by o. Then we proceed to compute explicitly the coupled equa¬ 
tions Xi = 0 and X 2 = 0, derive conserved quantities, and formulate 
a conjecture. The resulting equations are similar to those of Theorem 
(14.ip . except that m = —vn(xi — x 2 )(j)(t,xi)(j)(t,x 2 ) is replaced by 


0 = -v N (xi - x 2 ) 


^(t,x 1 )^(t,x 2 ) + — sh(2fc)(t,xi,a; 2 ) 


and similar O(jj) coupling corrections apply to the Hartree operator 
as well as S and W. 


Remark 5.2. The static terms of Xo(t) (not involving time derivatives) 
also appear in the recent preprint [T] , but do not serve as a Lagrangian 
there. 
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6. The Lagrangian and the equations, abstract 

FORMULATION 

Proposition 6.1. Let k and 0 be fixed. 


d_ 

de 


X 0 (<f) + eh, k) 

e=0 


(29) 


= 


Xi(t,x) ch (k)(t,x,y)h(t,y) - sh (k)(t, x,y)h(t,y)) dxdy 


In particular, if this vanishes for all h, then Xfit, x) = 0. 

Proof. 1-L re d can be written as 

nj _ 1 ^ p B p s/NA ( 1 ^ , nj \ _ — Vn A P —B 

n ’‘ i ~lWt +ee [~lm +n ) oe e 

in the sense of compositions (in space only) of operators. During this 
proof, denote hit = — + PL. 

Let h be an L 2 function and let 
A e = \f~N (a (0 + eh) — a* (0 + eh)). Thus we have 


A0(0 + eh, k) = 


e B e Ae 


H t e 


~Ae„~B I 


0>.|o> 


We compute 


d 

de 


e At ) e • 4(1 — Aq + — [Aq,Ao 


6=0 


d —Ao _ 


= V~N (a(h) — a*(h )) + — [a(0) — a*(0), a(h) 
= Vn (a(h) — a*(h )) + iiVA f <f>h 


a*(h)] 


e A 0 


d_ 

de 


e 


-Ae 


6=0 


d_ 

de 


e 


Ae 


6=0 


e~A° 


and 
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thus 


d 
de 

= { e 


/ e B e A t'H t e Af e B |o), |0> 

VNa(h) - VNa*(h),e Ao U t e- M le'^O), |0) 
e B (y/Na(h ) - VNa*(h)) e~ B , e 3 e^Ute^ e~ B |0>, |0) 


a(J) - a*(l), e B e Ao 'H t e~ Ao e~ B |0>, |o) 

= 2K^ red |0),a*(Z)|0>\ :=/ 
where we denoted 

e B (^VNa(h) — VNa*(h)^j e~ B = a(l) — a*(l) 

Explicitly, 

e B ( a(h ) — a*(h )) e~ B 
= a(ch(fc) o h) + a*(sh(fc) o h) 


a(sh(fc) o h) — a*(ch(k) o h) 


so 

l = V~N ^ch(fc) oh — sh (k) o hj 

Thus, 

/ =2V / iV3fj J Xi(t, x) ^ch(A;)(t, x, y)h(y) — sh(k)(t, x, y)h{y) s j dxdy 
=2 Vn$. J (ch(fc) o Xi - sh (k) o w) ( y)h(y)dy 

□ 


In order to state the corresponding result for X 2 , we have to intro¬ 
duce a new set of coordinates for our basic matrices 

k = ( M k ) sh(fc) \ 

\sh (k) ch(fc)/ 


where 

K __ ( 0 k(t , x, y)\ 

" \ k(t,x,y ) 0 J 


(30) 


The most obvious coordinate system is, of course, provided by k. We 
recall the following proposition, proved in [18], 
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Proposition 6.2. The exponential map is one-to-one and onto from 
matrices of the form (13Uj) (k 6 L 2 , symmetric) to positive definite 
matrices E satisfying the three properties of Remark (13. ip for which 
||/ — E \\ L 2 is finite. 


For our purposes, a better coordinate system is provided by ( = 
th(fc) = ch(fc) o sh (k). 

Proposition 6.3. There is a bijection between k e L 2 , symmetric, and 
( G L 2 , symmetric, ||C||o P < 1 (op stands for the operator norm) such 
that 


e 


k . 


= £<■■= 


f ch (k) sh(fc)\ 

\sh(fc) ch(fc )J 

(i c\(0-(°0 l/2 0_ 1 

i) \ 0 (/-C°CT 1/ V 

( (/-c°cr 1/2 co(/-coc)- 1/2 i 
bo(/-coo - 1/2 (/-c°cr 1/2 ) 




(31) 


where the square root is taken in the operator sense. 

Proof. Given k, define ( = ch (k) o sh(fc). The decomposition (j3B 
is an algebraic identity, and it is clear that ( is symmetric and L 2 . 
Since / — ch(/c)~ 2 = ( o we see that ||C||o P < 1- In fact, ||Cn||^ 2 = 
||n|ri 2 -||ch(fc)- 1 u|| 2 2 . Conversely, given ( asymmetric Hilbert-Schmidt 
kernel with ||C||o P < 1 define E^ by (l3B . It is easy to check that 
E^ is positive definite, satisfies the symmetries of remark (13.ip and 
||/ — E^\\hs < oo. (HS stands for the Hilbert-Schmidt norm), thus we 
can apply Proposition (16.21) and find the corresponding K. □ 


We also record the following consequence: 

Proposition 6.4. Let ( £ = ( + eh (h G L 2 , symmetric, ||C||op < 1/. 
and K f corresponding to ( e according to the previous proposition. Then 


d 

de 


e Ke e~ K 


e=0 



with 


b = ch [k) oho ch(/c) 
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Proof. We compute 
d 

e=0 


de 


e K *e~ K 


ch(fc) / o ch(fc) — sh(fc) o sh (k) — ch (k)' o sh(k) + sh(fc) o ch(fc) 
sh (k)' o ch (k) — ch (k) o sh (k) — sh (k) 1 o sh (k) + ch (k) o ch(fc) 


An easy calculation shows that b = —ch(fc) o sh (k) + sh(/c) / o ch(k) = 
ch(k) o ('o ch(k). □ 

We are ready to prove 

sx, 


— = —=ch(k) o X 2 o ch (k) 

Proposition 6.5. Let k e correspond to ( + eh as in the previous propo¬ 
sition. Then 


d 

de 


e=0 


X 0 (</>, kf) = a/ 25? / ch(fc) o X 2 o ch (k)(t, z, w)h{t, z, w)dzdw 


In particular, if the above vanishes for all h, then X 2 = 0. 
Proof. Let B e = B(k e ). 


_ / 0 B, 0 VnA tt c -VNA a -B, 


Xo(<j>, k e ) = (e^e' /NA H t e 


0>,|0> 


and 


d 

de 


( _ o A' 0 ( = -2s/i/« d |0),V>|0) 


(32) 


where 


ip = 


de 


e B *e- B = 1 


e=0 


d_ 

de 


K e —K 

e e 


e=0 


Using the isomorphism (J7|) and proposition (16.41) we see that 


^|0) = {id, 0,- —ch(k) oho ch.(k)(t, xi,x 2 ), 0, • • •) 

\/2 

where 6 is a real number coming from the trace of the self-adjoint 
a. Since X 0 is real, iO does not contribute to (1321) . and the result 
follows. □ 
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7. Explicit form of the Lagrangian 
The goal of this section is the following proposition. 

Proposition 7.1. The zeroth order term in TL re d |0) (which provides 
the Lagrangian density for our coupled equations) is X 0 (t) where 


X Q (t) = N / dx 1 \-$S (01 $0i) + | V 0! |" 


N 

~2 

1 

2 


+ / dxidx 2 v^_ 2 |0i02 + — (sh ° c h) 1,21 2 

dxidx 2 dx 3 Vi_ 2 \(j) 1 sh2 t 3 + 0 2 shi )3 | 2 


+ - ( / dxxdx 2 \ (shi )2 ^shi ;2 ) + |Vi, 2 sh lj2 | 


1 


2N 


dtdxidx 2 Vi_ 2 < |(sh o sh)i j2 p + (sh o sh)i i i(sh o sh) 2)2 


where sh 12 is an abbreviation for sh(k)(t, X\, x 2 ), etc, and the products 
are pointwise products, while compositions are denoted by o. 

The proof follows from several lemmas, which can be proved by ex¬ 
plicit calculations. We proceed to compute X 0 in (|25]1 . The only terms 
in (TTTll which contribute to X 0 are NVo which is already explicit, the 
zeroth order terms in T{R) |0), as well as the zeroth order terms in 
N-'ePVie- 3 10 ). 

Lemma 7.2. The term NVq is given by 


NV 0 = 


N 

N 

2 ~ 


dx \2 ^ ~ l V< ^ 

dxidx 2 {n^ 2 |0i0 2 | 2 } . 


We used abbreviations v(f 2 = vn(xi — x 2 ), 0i = 0(xi), etc., and 
for the following two lemmas we will denote u^ 2 = sh(k)(t, Xi, x 2 ) and 
ci, 2 = ch(k)(t,x 1 ,x 2 ). 
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Lemma 7.3. The zeroth order term inl(R) |0) is 

~\{j dx ^ dx ^ (“1,2^«1,2 - dtUi i2 u 1>2 ) + | Vi, 2 Ul, 2 | 2 

+ J dxxdx-idxz {u^ 2 |0iM 2 ,3| 2 + \<hvi, 3| 2 } 

+ 23? J dx ± dx 2 dx 3 

+ 23? / dxidx 2 {vi_ 2 (u ° c)i, 2 (j)i(j) 2 } 


Lemma 7.4. The zeroth order term in —j^e B Ve e |0) is 


2 / dx 1 dx 2 v^_ 2 \ (tt°c)i i2 (ttoc)i, 2 


+ |(ttO It) l t 2 | 2 + (u o lt)l j l(lt o u) 2 ,2 


8. Explicit form of the equations 

In this section we derive the following theorem, thus introducing our 
new equations. First, some notation. Consider the kernels 

w c (t,x,y) =j>(t,x)(f>(t,y) 

u p (t, X, y) = sh (k) O sh (k) (t, x, y ) 

and their trace densities 

Pc = \<P\ 2 (t,x ) 

p p (t,x ) = sh (k) o sh (k)(t,x,x) 

Here c stands for condensate, and p for pair. In this notation, the old 
operator kernel defined in (j20c[) is 

g N (t, x, y) := -A x S(x-y) + (v N * p c )(t,x)S(x - y) 

+ v N (x~y)u c (t,x,y) 

Dehne the new operator kernel 

g N {t,x,y) := -A J(x - y) 

+ (v N *p c )(t,x)S(x-y) + v N (x - y)u c (t , x, y) (33) 

1 

+ (K * Pp){t, x)S(x -y) + v N (x - y)u p (t , x, y)) (34) 
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Also denote a c = (j33j) , j^a p = (1511) and a = a c + j^a p . Define 
S(s) := -s t + g%os + sog N and W(p) : = -p t + [gj,,p\ 

% i 

Finally, define @(f, xi,x 2 ) = — njv(t, Xi, x 2 ) (4>(t, Xi)<f(t, x 2 ) + j^sh(2k){t, Xi, x 2 )) ■ 
Theorem 8.1. The equation Ad = 0 is equivalent to 

Q(t,x 1 ,x 2 )(j)(t,x2)dx2 + j -^a p (t,x 1 ,x 2 )(j)(t,X2)dx2 = 0 

The equation X 2 = 0 is equivalent to either of : 

1) the equation 

S(th(/c)) = 0 + th(fc) o 0 o th(fc) 

2) the pair of equations (in fact, 2a) implies 2b)) 

2a) S (sh(2/c)) = 0 o ch(2fc) + ch(2 k) o 0 (35) 

2b) W (ch(2kj) = © ° sh(2fc) - sh(2k) o 0 

Remark 8.2. One can go back and fourth between ( and ch(2fc), sh(2fc) 
using 

sh (k) o sh (k) — (l — Co C) _1 — 1 = - (ch(2fc) — 1) 

( = sli(2/c)(l + ch(2/c)) _1 
Proof. A direct calculation for X { shows that 

Ad = —VN ^ch(fc) o Har k ((j)) + sh {k) o Har k ((J. >) j 

where 


^d t (j)(t,x i) - Acf) - J 


Har k (<t>)(t,x i) 


= -d t 4> - A0 

i 


Q(t, xi, x 2 )4>{t, x 2 )dx 2 


+ J VN^Xx - x 2 )(shosh.)(x 1 ,x 2 )(l)(x 2 )dx 2 
+ jj<t>{x i) J v N (xi — x 2 )(sh o sh)(a: 2 , x 2 )dx 2 

In conjunction with Proposition (16. ip this shows that 

= NHar k {4>) 


which can also be easily verified directly from Proposition (17.11) . 
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A direct calculation also shows that, if X 2 denotes the second com¬ 
ponent of T-Lredl 0), then 


V2X,(t,y h y 2 ) = 


(36) 


S(sh(k)) — ch(k) o raj o ch(A;) — (W(ch(£;)) + sh(A;) o mj o sh(A;) 


+(1/N) / dx\dx 2 


ch(yi, x 2 )sh(x 2 , 1 / 2 ) (sh o sh) (aq, x\)vn(xi - x 2 )+ 

ch(yi, x 2 )sh(aq, y 2 ) (sh o sh) (aq, x 2 )v N (xi - x 2 )+ 
ch(j/i, xi)sh(x 2 , y 2 ) (sh o sh) (aq, x 2 )v N (x 1 - x 2 )+ 

ch(yi, aq)sh(aq, y 2 ) (sh o sh) (x 2 , aq)v/v(aq - x 2 )\ + 

/ symm 

sh(j/i, aq)sh(aq, y 2 ) (sh o ch) (aq, x 2 )v N (x 1 - x 2 )+ 
ch(t/i, aq) ch(x 2 , y 2 ) (ch o sh)(aq, x 2 )vn(x\ — x 2 )| . 

where symm stands for ” symmetrized”. The time dependance in the 
last six lines has been omitted. Recalling ( = ch(fc) osh(k) = sh(k) o 

ch(k) , compose on the left with ch(k) and on the right with ch(/c)^ 1 

to get 


—1 


1 


ch(Jfe) oXjoch(fc)- 1 = S(C)-©-C°0°C + ^N (37) 


where N is given by 


N(t,yi,y 2 ) =C(t,yi,y 2 )( / dx ( (sh o sh + sh o sh) (t, x, x)v N (x± Vl ) J + 

\ J \ / symm 

f dx((t, x, y 2 ) (sh o sh + sh o sh) (t, x, i/i)ujv(x — y\) 


symm 


where symm stands for symmetrizing in y\ , y 2 . In other words, 

N = (oa p + tt p o( 

Thus, in ( coordinates, the equation X 2 = 0 becomes 


s(o-e-c°e°C = o 


( 38 ) 
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Now we can get an equation for W(ch(2/c)) and S(sh(2 k)). We will use 
the general formulas 

W(/- 1 ) = -/- 1 oW(/)o/- 1 

W(f°g) = S(f)og-foS(i) 

s (/ o #) = S (/) O g - f oW (g) 


Thus 


w((i-CoC)- 1 ) = (1 -CoC)- 1 o(s(0oC-C°s(C)) 

= (i -C oC)- 1 °((0 + C°0°C)C + C°(S + C°0°C)')(i -c°o _1 


Similarly we get a formula for S(sh(2/c)), using 

= (i-C°C)" 10 (s(C)-C°S(C)°c) °(i-C°C ) _1 


s(C°(i-C°C) *) = 


(i-C°C) 1q @ + C 0 ©°C + C° (© + C°@°C)°C °(i~C°C) 1 


= ( 1 -C°C) - 7T 00 + 00 (l-(o() - 


□ 
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9. Conserved quantities 

We start by motivating the introduction of some conserved quanti¬ 
ties. Recall the Lagrangian 


C((j),sh(k)) — N j dtdxi |—^ (4>id t (j)i) + | V0i | 

N f 1 

+ — dtdxidx 2 Vi_ 2 \(j) 1 (f >2 + — (sh o ch)i i2 | 2 

+ ^ J dtdxidx 2 dx 3 Vi_ 2 \(/) 1 sh.2,3 + 0 2 shi, 3 | 2 

( j dtdx l dx 2 | (sh 1)2 ^shi i2 ) + |Vi i2 shi j2 | 2 


+-— / dtdxidx 2 Vi_ 2 ^\(sh o sh)i j2 | 2 + (sh o sh)i j i(sh osh) 2)2 


where sh 12 is an abbreviation for sh (k)(t, xi, x 2 ), etc, and the products 
are pointwise products, while compositions are denoted by o. Introduce 
the energy 8 


8(<f>,sh.(k))(t) —N dx i I |V0i 


+ y I dxidx 2 v^_ 2 \4>\4> 2 + ^(sh o ch)! 2 1 2 

+ \ [ dxidx 2 dx 3 Vi_ 2 \(f)ish. 2t3 + 0 2 shi ;3 | 2 


+-( / dx\dx 2 || Vi i2 sh 12 | 2 | 


+ -— / dxidx 2 v^_ 2 )y |(sh o sh)i i2 | 2 + (sh o sh) L1 (sh o sh) 2>2 


Our equations for 0 and sh(k) are equivalent to 


1 <90 88 

i d t 8~4> 

1 dsh(k) _ 88 

* dt 8sh(k) 


( 40 ) 
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The relation 


d 


0 =^\ e=0 £(e l <t>,e 1 sh (k)) 
=2!R [ / ^=(— i(f))dx\ + j 


58 


—ish(k))dxidx2 I 

5sh(k )' J 

together with (j39jl . (l40lh leads to the conservation 

— ( [ |0(t, xi)\ 2 dx\ + — [ |sh(£;)(t, x\, X 2 ) \ 2 dx\dx 2 J =0 


dt 


thus we define the density 


N 


p(t,X!) = | 0 (t,Xi )| 2 + — J \sh(k)(t,x 1 ,x 2 )\ 2 dx 2 
= p c (t,x 1) + ^p p (t,x 1) 

Similarly, let 4> e (t,x) = (j)(t, x + ee 3 ), sh(fc) e (t, x, y) = sh(fc)(t,x + 
eej, y + eej) (ej = unit vector, 1 < j < 3). The relation 


d 

°~Te 


e=0 


8(cj) e ,sh(k) e ) 


SS 


=23? ( / j=d 3 (pdx\ 


58 


5sh(k) 

together with (13U1) . (jlOj) leads to the conservation 


( djsh(k))dxidx2 


d_ 

dt 


N J ^ [4>dj4>) dx 1 + J ^ ^sh(fc)9jsh(fc) j dx\dx^j = 0 


thus we define the momentum density 

1 


Pj(t,x 1 ) =9‘ ^sh(fc)9jsh(A;)^ dx 2 


Finally, using 


:= Pc,j(t, x i) + — 


we see that 


= 0 
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so we define the energy density 
e(t,x i) = iV|V</>i|“ 

+ y I dx 2 Vi_ 2 \(j> 1 (i>2 + -^(sho ch )i,2| 2 

+ ^ J dx 2 dx 3 Vi_ 2 10ish 2 , 3 + 02shi j3 | 2 
J dx 2 11Vi,2shi j2 1 | 

+^/^<- 2 {| ( shoih )li2 r + ( S hoih )l yihosh k2 }) . 

10. A CONJECTURE 

We conjecture that, if </>, fc satisfy the equations of Theorem (18.11) 
and | ^exact), | 'ipappr )j are defined by (HU), (TT5D . then, in the critical case 
/3 = 1, 

|||f/Yrarf) |l/7ippr^ ||T ^ 0 

as N —y oo, at an explicit rate. 
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